Stability of the doped antiferromagnetic state of the t — t'-Hubbard model 



Avinash Sine and Haranath GhosbEl 
Department of Physics, Indian Institute of Technology Kanpur - 208016, India 

The next-nearest-neighbour hopping term t' is shown to stabilize the AF state of the doped Hub- 
bard model with respect to transverse perturbations in the order parameter by strongly suppressing 
the intraband particle-hole processes. For a fixed sign of t', this stabilization is found to be signif- 
icantly different for electron and hole doping, which qualitatively explains the observed difference 
in the degree of robustness of the AF state in the electron-doped (Nd2- x Ce x Cu04) and hole-doped 
(La2- x Sr x Cu04) cuprates. The t' — U phase diagram is obtained for both signs of the t' term, 
showing the different regions of stability and instability of the doped antiferromagnet. Doping is 
shown to suppress the i'-induced frustration due to the competing interaction J'. A study of trans- 
verse spin fluctuations in the metallic AF state reveals that the decay of magnons into particle-hole 
excitations yields an interesting low-energy result F ~ ui for magnon damping. 



I. INTRODUCTION 

The problem of determining the self-consistent mean- 
field (MF) state of an interacting electron system is 
generally equivalent to minimizing the ground-state en- 
ergy function Eg{A} with respect to an assumed order- 
parameter set {A}. However, it may happen that when 
this set is expanded to include an additional component, 
the original self-consistent fixed point turns out to be 
a local maximum with respect to this new direction in 
order-parameter space. In other words, the original self- 
consistent fixed point actually represents, in such a case, 
a saddle point in the expanded order-parameter space. 

This is precisely the fate of the antiferromagnetic 
(AF) state of the doped Hubbard model with nearest- 
neighbour (NN) hopping t. The self-consistent AF state 
is described by the gap parameter A z , related to the 
sublattice magnetization m z in the z direction by 2A Z — 
m z U . This state was found to be unstable with respect 
to a new direction represented by the transverse compo- 
nent Aj_. Inclusion of the transverse magnetization m± 
amounts to spin twisting away from the perfect AF align- 
ment, and the above instability is, in fact, towards the 
so-called spiral state, [jO-l!) characterized by an ordering 
wave vector Q, different from the AF wave vector (tt, 7r). 
The spiral state is stabilized by an energy gain (for the 
doped system) resulting from an activation of the 0(t) 
hopping due to the slight twisting of neighbouring spins 
relative to the AF alignment. In terms of the devia- 
tion q = Q — (tt, 7r) from the AF wave vector, the twist- 
induced hopping energy gain ~ tq was found to compete 
with the twist-induced exchange energy loss ~ Jq 2 , re- 
sulting in an optimum spiral pitch for a given doping 
concentration. Thus, minimization of Eg{A}, after ex- 
panding the order-parameter set to include the transverse 
component Aj_ (or equivalently the spiral wave vector 
Q), results in a fundamentally new self-consistent fixed 
point. 

In this paper we show that for a positive next-nearest- 
neighbour (NNN) hopping term t' in the Hubbard model 
(see Eq. 1 for our convention), the above instabil- 
ity is strongly suppressed for hole doping, and the self- 



consistent AF state of the doped t — t'-Hubbard model 
becomes stable for a range of doping concentration. || 
However, for negative t' and hole doping, we find that 
while the AF state survives longitudinal perturbations 
in the order parameter, it remains unstable with re- 
spect to transverse perturbations for any doping. To 
the extent the t — i'-Hubbard model is equivalent to the 
three-band Cu-0 Hamiltonian, this result is par- 

ticularly relevant to doped cuprates in which AF order- 
ing is much more robust with respect to electron dop- 
ing. The electron doped cuprate Nd2- x Ce x Cu04 retains 
AF order up to a doping concentration of about 12%, 
B whereas only 2% hole concentration destroys AF or- 
der in La2- x Sr x Cu04. Since a particle- hole transforma- 
tion maps the t' model and hole (electron) doping on the 
—t' model and electron (hole) doping, the positive (nega- 
tive) t' model and hole (electron) doping is appropriate to 
study for the electron-doped compound Nd2- x Ce x Cu04, 
for which a negative t' is usually assumed. 

More significantly, a stable AF state of the doped t—t'- 
Hubbard model provides a microscopic realization of the 
metallic AF state, in which the Fermi energy lies within a 
quasiparticle band. Indeed, metallic antiferromagnetism 
has been reported in k — (BEDT — TTF) 2 X, § V 2 - x 3 , 
JlO| and NiS2- x Se x . jllj This allows us to quantitatively 
study spin excitations in the doped AF. Specifically, we 
evaluate the magnon damping T arising from magnon 
decay into intraband particle-hole excitations across the 
Fermi energy. Interestingly, we find that in the low- 
energy limit, the magnon damping r ~ u>. This is of rele- 
vance in the context of phcnomenological theories such as 
the nearly AF Fermi-liquid theory, Jl2[ put forward to ex- 
plain the Knight shift experiments and describe the non- 
Korringa temperature dependence of the nuclear spin- 
lattice relaxation rate in doped cuprates. 

The need for more realistic microscopic models which 
include NNN hopping etc., has been acknowledged re- 
cently from band structure studies, photoemission data 
and neutron-scattering measurements of high-T c and re- 
lated materials. Jil^dTq] Estimates for \t'/t\ range from 
0.15 to 0.5. Effect of hole and electron doping on the 
commensurate spin ordering have been studied for the 
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t — i'-Hubbard model and applied to La2- x Sr x Cu04 and 
Nd2-xCe x CuC>4. 17 Spin correlation function, incom- 
mensurability, and local magnetic moments in the doped 
t — i'-Hubbard model have been studied using the Quan- 
tum Monte Carlo method. jl8| At half filling, existence of 
a metallic AF phase has been suggested in d — 2. jl9| The 
suppression of the perfect-nesting instability by the NNN 
hopping, and the critical interaction U c vs. t' phase dia- 
gram has been studied in d — 2, 3. |2(]] Magnon softening 
due to t' and a significant enhancement in the low-energy 
spectral function due to single-particle excitations has 
been observed. ]2l[ | Evolution of the magnon spectrum 
with if was studied, and the quantum spin-fluctuation 
correction to sublattice magnetization in d — 2 and the 
Neel temperature in d = 3 were also evaluated. J2l]] 
Finite-?/ effects on competing interactions and frustra- 
tion were recently examined, and the magnetic phase di- 
agram was obtained at half filling. [^2| Recently effects of 
nnn hopping have also been studied in one-dimensional 
systems involving chains and ladders. fl23[] 



II. METALLIC AF STATE 

We consider the t — i'-Hubbard model on a square lat- 
tice, with NN and NNN hopping terms t and t' connect- 
ing sites i to i + 5 and i + k, respectively: 
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In the following we set t = 1. As the t and t' terms 
connect sites of opposite and same sublattice, respec- 
tively, the corresponding free-fermion energies e k = 
— 2t(cosk x + cos k y ) and e k = —4t' cosk x cosk y appear 
in off-diagonal and diagonal matrix elements of the HF 
Hamiltonian, and in the two-sublattice basis we have 
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for spin a. Here 2A = mil, where m is the sublattice 
magnetization. Since the e k term appears as a unit ma- 
trix, the eigenvectors of the HF Hamiltonian remain un- 
changed from the NN hopping case. The fermionic quasi- 
particle amplitudes on the two sublattices A and B are 
given by the eigenvector (a k(J 6 kcr ) of the HF Hamilto- 
nian matrix, and for spin a =t, J, and the two quasipar- 
ticle bands 0, ®, we have IM 
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FIG. 1. The AF-state electronic density of states (DOS) 
for positive t 1 (solid) and negative t' (dashed), showing the 
flat DOS near the Fermi energy for hole and electron doping, 
respectively. 

In the strong coupling limit (U >> t), the majority and 
minority densities given above reduce to 1 — e£/4A 2 and 
e k /4A 2 , respectively. 

The quasiparticle energies corresponding to the eigen- 
values of the HF Hamiltonian are, however, modified to 
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the two signs ± referring to the lower and upper quasipar- 
ticle bands. The band gap is thus affected by the NNN 
hopping term, and it progressively decreases as 2A — 4t' 
in the weak coupling limit. The corresponding density of 
states is shown in Fig. 1 for both positive and negative t', 
showing the drastic reduction near the upper and lower 
band edges of the lower and upper bands, respectively. 

We first consider the case of positive t' and hole dop- 
ing; the same results hold for negative t' and electron 
doping. The highest electronic energy levels of the 



\J A 2 + e 2 . correspond to 
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lower AF band 

(k x ,k y ) — (±7r, 0) and (0,±7r), and with doping hole 
pockets develop around these points. In this region of 
k space constant-energy surfaces are nearly (semi) circu- 
lar, as shown in Fig. 2. Therefore, for low hole doping 
the Fermi surface consists of these four semi circles. The 
radius a of the hole Fermi circle is related to the dop- 
ing concentration x. The ratio of the total area inside 
the Fermi circles (27ra 2 ) to the total Brillouin zone area 
(47r 2 ) yields the fraction of unoccupied states (x), so that 
a 2 = 2nx. 

We now examine the AF states which lie near the Fermi 
energy E-p, and will therefore contribute to the particle- 
hole processes at low energies. For k states lying close to 
a Fermi circle, say near the point (— 7r, 0), k Xl k y can be 
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2. Constant energy surfaces for the AF band energy 
k — (A 2 + tk) 1 ^ 2 , showing the nearly (semi) circular 
contours near the top of the lower energy band around the k 
points (±7T,0) and (0, ±tt). Here t' = 0.05 and A = 5. 
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parametrized in terms of the angle 8 shown in Fig. 3, 
so that k x = — (tt — a cos 9) and k y = asin#. Thus, 
cos k x ps -1 + (a 2 /2) cos 2 6* and cosfcy pa 1 — (a 2 /2) sin 2 
up to 0(a 2 ), so that ek = — 2i(a 2 /2) cos 26* ~ tx whereas 
e' k = 4f'(l — a 2 /2) ~ t'. Therefore, in the strong coupling 
limit the AF band energy E^ pa e k — e k /2A — A, which 
further simplifies to e k — A in the limit Jx 2 « t'. This 
yields a nearly quadratic energy dispersion for k close to 
(±7r, 0) or (0, ±7r), which accounts for the nearly circu- 
lar constant-energy surfaces and the flat density of states 
near the top of the lower band, as seen in Fig. 1. 

For negative t' and hole doping, the hole pockets in 
the lower band develop around the points (±7r/2, ±7r/2). 
Translating the k-space origin to (tt/2,tt/2), and then 
rotating the coordinate system by 7r/4, the quasiparticle 
energy E^ can be written in terms of the transformed 
coordinates k as 
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to second order in k x ,k v , indicating that for J > \t'\ 
the constant-energy surfaces are ellipses centered at 
(7r/2,7r/2). Here the strong coupling limit has been 
taken, and the band energy shifted by A. Equating the 
total normalized area of the four elliptical hole pockets 
to the doping concentration x, we obtain the semi- major 
and semi-minor axes of the Fermi-ellipse 
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yielding the same y/x dependence as obtained earlier for 
the Fermi-circle radius for positive t' . 



III. STABILITY OF THE AF STATE 



To test the stability of the AF state with respect 
to transverse perturbations Aj_ in the order parame- 
ter, we examine the transverse response within first- 
order perturbation theory. For a transverse perturbation 



(n,0) {A^}, which leads to terms A^_aLa^ + (A^)*a^a 4 | in 
the Hamiltonian, the first-order corrections to the wave- 
functions of the HF Hamiltonian H(A Z ) are first ob- 
tained. Next, the transverse mean-field potential Aj_ = 
— Z7(at, Oif) developed due to these corrections is ob- 
tained upto first order in A^; this may be expressed 
as Aj_ = UJ2j[^]ij^±^ where the transverse response 
matrix [A] connects the cause and effect. If the largest 
eigenvalue A™ ax of this matrix [A] is larger (smaller) than 
l/U, then within a self-consistent approach, the AF state 
characterized by the single order parameter A 2 is unsta- 
ble (stable) with respect to spin twisting. 



It is shown in the Appendix that the transverse re- 
sponse matrix [A] is nothing but the bare antiparallcl- 
spin particle-hole propagator [x°] in the static limit (lu = 
0). In the AF state, where translational symmetry holds 
within the two-sublattice basis, Fourier transformation 
to wave vector (q) space reduces [x°] to a 2 x 2 matrix 
[x (q)] having two eigenvalue branches, which are further 
discussed below. 

For the undoped t — i'-Hubbard model in the strong 
coupling limit, the AF insulating state involves a com- 
pletely filled lower band and empty upper band, so 
that only interband particle-hole processes contribute 
to [x°(q)]- The relevant eigenvalue branch (having 
the larger eigenvalues) is given by A 1 f ltcr (q) = 1/C — 
(t 2 /4A 3 )(l-2J'/J)g 2 for small q. || The corresponding 
eigenvector is (1 — 1), indicating that the transverse fluc- 
tuations AY = (— I) i A'^e lqr for the mode q involve op- 
posite sign on the two sublattices. This implies that the 
mode q simply represents a spiral twisting of spins with 
an appropriate wavelength. For J' < J/2, the eigenvalue 
is less than l/U for all finite q, indicating the stability 
of the AF state, whereas for J' > J/2, the AF state is 
unstable towards a F-AF state with ordering wavevector 
Q = (7t, 0), possessing antiferromagnetic (ferromagnetic) 



ordering in the x (y) direction. 1 22 1 The q = mode corre- 
sponds to a rigid rotation of al 
response for this mode (C/Aj 



1 spins, and the marginal 
1) merely reflects the 
equivalence of all spontaneous symmetry-breaking direc- 
tions, which follows from the Goldstone theorem for sys- 
tems possessing continuous symmetry. 

With doping the Fermi energy moves into a quasipar- 
ticle band, and particle-hole excitations now include in- 
traband processes involving states of the same band. By 
including these intraband processes in the transverse re- 
sponse matrix [x°], it was shown earlier for the doped 
Hubbard model that the eigenvalue increases beyond l/U 



3 



for finite q, indicating that the AF state is unstable to- 
wards the spiral twisting of spins for any finite doping. (^] 
As the Hubbard model exhibits particle-hole symmetry, 
this instability applies to both hole and electron doping. 

In this section we show that in the t — i'-Hubbard 
model this intraband contribution is strongly suppressed 
for positive (negative) t' and hole (electron) doping, re- 
sulting in a stabilization of the AF state for a range of 
doping concentration. This suppression is due to the 
modification in the quasiparticle energy spectrum by t', 
which leads to a significant reduction in the phase space 
available for intraband processes. However, for negative 
t' and hole doping, the AF state is found to remain unsta- 
ble with respect to transverse perturbations in the order 
parameter for any doping. 

In order to obtain the transverse eigenvalue Aj_, 
we evaluate, in the static limit (u> = 0), the bare 
antiparallel-spin particle- hole propagator, [x (q, ui)] = 
i/^£ k [G T (W)][G±(k-q,a/ - u)] in the AF state, 
separately focussing on the intraband and interband pro- 
cesses. For the case of hole (electron) doping, intraband 
processes involve particle and hole states from the lower 
(upper) band, while interband processes involve particle 
and hole states from different bands. 

Now, the matrix [x°(q)] has two eigenvalue branches 
[A±(q)] 1(2) = [x°(q)]AA ± [x°(q)]AB, corresponding to 
the two eigenvectors (1 1) and (1 — 1), respec- 
tively. Separating the interband and intraband contri- 
butions as [x°(q)] = [x°(q)] intor + [X°(q)] intra , the total 
transverse eigenvalues can be written as [A^(q)] 1 ^ 2 -' = 
[A i f cr (q)] 1(2) + [A^f tra (qi)] 1(2) . We focus on the (1 - 1) 
branch which yields larger eigenvalues for small q, even 
when intraband processes are included. 

In the following we show that for the (1 — 1) 
branch the intraband contribution can be written as 
A^ tra (q) = (i 2 /4A 3 )a intr V for small q, where a intra 
is a dimensionless positive coefficient. The interband 
contribution can similarly be written as A^ ter (q) = 
l/U — (t 2 /4A 3 )a mter 5 2 . The total transverse eigenvalue 
is therefore obtained as 



A 1 (q) = Ar(q)+Ar ra (q) 
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For q = the transverse eigenvalue exactly equals l/U, as 
required from the Goldstone theorem and the continuous 
spin-rotation symmetry. Thus the question of stability 
of the AF state with respect to transverse perturbations 
is reduced to a comparison of the two coefficients a mtra 
and a mter . If a lntra exceeds a mter , the transverse response 
UX±(q) > 1 for finite q, signalling a spiral instability. 



A. Intraband contribution 

In the two-sublattice basis (labelled by A,B), the in- 
traband contributions to the diagonal and off-diagonal 



elements are obtained as 
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where £ k indicates that states k, — k are below the Fermi 
energy E-p, while states k — q, — k + q are above E-p. 

Simplification results from the property that the band 
energies -E k and the amplitudes (a k b k ) are both un- 
changed under the transformation k — > — k. The strong 
coupling limit U » t leads to further simplification 
and in the static limit (ui — 0) we obtain for the diag- 
onal [x°(q)]d = [x°(q)]AA = [x°(q)]BB and off-diagonal 
[x°(q)]od = [x°(q)]AB = [x°(q)]BA matrix elements 



[x Q (q)] d ntra = E 77 
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k (4- q -^ k )-(4- q -4)/2A 
[x°(q)]od tra = E 7? iTSui^ 3^T - 
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From Eq. (9) the intraband eigenvalue A™ tra (q) = 
[x°(q)]AA a - [x°(q)]AB a for thc relevant (1 - 1) branch 
is obtained as 



*(q) = E 
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The relabelling transformation k x — > — k x and/or k y — > 
—k y must leave A™ tra (q) invariant, and as both e k and 
e k are even functions of k x and k y , A™ tra (q) must also 
be an even function of q x and q y . Furthermore, since the 
x and y directions are equivalent, A™ tra (q) must vary as 
q x + q y to lowest order in q. As Aintra (q) depends only 
on the magnitude q, we consider q y — and q x = q for 
simplicity. Figure 3 shows the hole Fermi circle of radius 
a near the {—it, 0) point in the Brillouin zone. Consider 
another semi-circle whose center is shifted by the wave 
vector q, as shown. For k states lying in the crescent 
region between these two semi-circles, we have E^ < Ep 
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FIG. 3. The k states which contribute to the intraband 
particle-hole processes near the (— 7T, 0) point lie in the cres- 
cent-shaped region. 

and -E^-q > Ep , and therefore it is this region of k space 
which contributes to the intraband processes in Eq. (10). 
As discussed earlier, for q << a, the Fermi circle radius, 
these k states can be parametrized as k x = — (n — a cos 9) 
and k v = asm.9. As a 2 = 2irx, for small doping concen- 
tration we have a << 1. 

We now show that the two terms in the energy de- 
nominator in Eq. (10) arc in the ratio t' : Jx, which 
plays a key role in determining the intraband contribu- 
tion. Taylor expansion to first order yields £k- q — £k = 
— q.Vek = — 2tqsink x and e k _ q — e k = — q.Ve k = 
— 4t'q sin k x cos k y for small q. Therefore, we have £k- q — 
£k = 2tqsn\(acos9) w 2tqacos9 for a << 1 and 
e k- q — e k = 4i'gsin(acos#) cos(asin#) « At'qacosO. As 
£k ~ tx, it follows that the ratio of the two terms in the 
energy denominator is ~ t' / Jx. In the t' « Jx limit, 
when the NNN hopping term may be dropped in com- 
parison, it was shown that A™ tra (q) actually diverges as 
x — ► 0, so that the AF state is unstable for any doping 
concentration. p| 

However, for t' » Jx the exchange-energy term may 
be dropped in Eq. (10). Integrating over the crescent- 
shaped region shown in Fig. 3, and multiplying by 2 to 
account for the additional contributions from the regions 
near (0,±7r), we obtain 
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Comparison with Eq. (7) yields cv mtra = Ax/2t', which 
vanishes as x — * 0, in sharp contrast to the result for 
t' « Jx. A similar result for a mtra proportional to x 
was obtained earlier. Q 



To evaluate A™ tra (q) for negative t' and hole doping, 
we take q x = q y for simplicity, so that q = (V2q, 0) in 
the momentum coordinate system rotated by 7r/4. Sub- 
stituting - E§ w 4V2(J - \t'\)K x q from Eq. (5) 
and £k- q — £k ~ 4tq in Eq. (10), we obtain 
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As the phase space available for the intraband process is 
of order ag, and k x is of order b, we obtain 



, intra 



V(J-W\)\t'\ 



(13) 



In the strong coupling limit a mtra >> 1, and is therefore 
much larger than cv lnter ~ 1, so that the transverse re- 
sponse eigenvalue increases with q, and the AF state is 
unstable for any doping concentration, as also found in 
ref. [17]. This result is confirmed in the numerical study 
described in section IV, where it is seen that in the t' — U 
region where the AF state exhibits longitudinal stabil- 
ity, the maximum transverse response eigenvalue remains 
well above 1/U for all U. 



B. Interband contribution 

The interband eigenvalue A™ ter (q) is next obtained 
for t' > by evaluating the interband contribution to 
[X°(q)], involving k states such that < Ep and 



> Ep. The second condition holds for all k in the 
strong coupling limit, and the first condition is met by 
excluding in the k sum states lying within the four Fermi 
circles near (±7r, 0) and (0, ±7r), as shown in Fig. 2. Fol- 
lowing the evaluation of the matrix elements of [x°( c l)] f° r 
the undoped t — t'-Hubbard model in the strong coupling 
limit, p2] | and evaluating the contribution from within 
the Fermi circles (radius a) up to order a 4 , we obtain 
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The off-diagonal elements are unchanged to order a 4 , the 
leading correction being of order a 6 . 
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Retaining terms only up to 0(a 2 ) and 0(q 2 ), we obtain 



A^ tGr (q) corresponding to the (1 
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where the doping concentration x = a 2 /2ir. In the un- 
doped case, the q 2 term changes sign when J' > J/2, 
reflecting the instability of the AF state towards the F- 
AF state with Q = (ir, 0). Doping is seen to effectively 
reduce this frustrating effect of t'. 

The phase boundary of the region of stability of the AF 
state, obtained by equating the two coefficients a mtor and 
tfntra from Eqs. (11) and (15), is therefore given by 
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FIG. 4. Plot of the function f(t') = (At' 3 - 2t')/(At' 2 - 1) 
defined in Eq. (17) showing the two branches. Intersec- 
tions with Ax yield the two bounds 1? min and ^axi between 
which the AF state is stable. For Aa; = 0.8, t^m ^ 0.3 and 
/ -11 

U 111 flX ~ J-.-L. 
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This is a cubic equation in t' which may be written as j-y. PHASE DIAGRAM 



, „ At' 3 - 2t' 

m = = . (l?) 

The function f(t') has two branches, as shown in Fig. 4, 
yielding two solutions for any Ax. The first branch in- 
creases as 2t' for t' << 1, and diverges as t' — > 1/2. The 
second branch starts from t' = 1/ \/2, and asymptotically 
approaches i' for t' » 1. This implies that for small Ax, 
the two solutions are i^ in w Ax/2 and t^ ax w 1/V2, and 
the AF state is stable for t' min < t' < t' max . Furthermore, 
t' min is bounded by 1/2 from above, t' max is bounded by 
1/V2 from below, and t' max asymptotically approaches 
Ax for large Ax. 

To summarize, for fixed x, the minimum il required 
to stabilize the AF state increases in proportion to the 
interaction strength as t' min w Ux/4, and for fixed t', the 
instability to a transverse perturbation develops when 
x > x c , where the critical doping concentration x c = 
2t'/A ps At' /U. Moreover, as t' max increases away from 
l/v2 with increasing doping, we find that hole doping 
interestingly suppresses the i'-induced frustration and 
stabilizes the AF state against the instability towards 
the F-AF state. These qualitative behaviours are indeed 
confirmed in the numerical study described in the next 
section. For a quantitative comparison we take a spe- 
cific case x = 0.09 (18 holes in a 14 x 14 lattice), and 
A = mU/2 a (1 - a) £7/2 « 0.9 for £7 = 20, so that 
Ax w 0.8. Graphical solution of Eq. (17) by consid- 
ering the intersections in Fig. 4 yields t' min ~ 0.3 and 
'•max ~ 1-h m excellent agreement with the numerical 
result. 



A collinear mean-field state such as the AF state, with 
spin ordering along some fixed (say, z) direction, can be 
unstable with respect to either longitudinal or transverse 
perturbations. Longitudinal perturbations modify the lo- 
cal magnetization (S l z ) and therefore the local gap A*, 
whereas transverse perturbations induce spin twisting, 
resulting in transverse terms A^_. In this section we de- 
scribe a numerical procedure for studying the stability of 
the AF phase with respect to both these perturbations, 
and discuss results for a 14 x 14 lattice. The objective 
is to obtain the t' — U phase diagram for a fixed dop- 
ing concentration, showing the different regions of stabil- 
ity/instability of the doped antiferromagnet. 

The usual numerical self-consistency procedure (with 
a spin-diagonal mean-field Hamiltonian corresponding to 
ordering in the z-direction) is implicitly sensitive to insta- 
bility with respect to longitudinal perturbations due to 
the presence of small noise associated with the numerics. 
Thus stability with respect to longitudinal fluctuations is 
automatically indicated, at least for the finite-size lattice, 
if a homogeneous AF state is obtained self-consistcntly. 
Instability with respect to transverse perturbations is 
studied by numerically evaluating the largest eigenvalue 
A" lax of the transverse response matrix [x°] ■ Numerical 
evaluation of [x ], using the electronic eigenvalues and 
eigenf unctions of the self-consistent AF state, has been 
described earlier. [Q The largest eigenvalue A™ ax exceed- 
ing 1/E7 signals the transverse instability. 
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FIG. 5. The t' — U phase diagram for positive t' and hole 
doping of around 9%, showing the region of stability between 
t' min and imax, where the AF state is stable with respect to 
both longitudinal (||) and transverse (±) perturbations. 



A. t' > 

Figure 5 shows the t' — U phase diagram for a 14 x 14 
lattice with a fixed doping concentration of around 9% 
(nine holes per spin). The minimum value t' min above 
which the AF state is stabilized is indeed found to in- 
crease nearly linearly with U in the strong coupling limit, 
and asymptotically approaches 0.5 with increasing U. 
The stable region (with respect to both perturbations) 
lies in a band between t' min and t' max . The instability 
for t' > t' maK is towards the Q = (it, 0) state, and t' max 
increases away from l/y/2 in the strong coupling limit. 
These features are in agreement with the analytical re- 
sults of the previous section. With further increase in t' 
the band gap continually decreases, causing an enhance- 
ment in the longitudinal response. This eventually re- 
sults in a longitudinal instability at t' — t[ ong where the 
gap vanishes. As expected, t[ ong scales linearly with U in 
the strong coupling limit. Longitudinal instability is also 
seen for t' ^ 0.03, essentially independent of U. 



B. t' < 

For negative t' the AF-state energy eigenvalues in the 
lower band are four-fold degenerate, therefore the num- 
ber of doped holes is taken in multiples of four. The 
region of stability of the AF state with respect to lon- 
gitudinal perturbations is shown in Fig. 6 for 8 doped 
holes per spin (x w 8%). For t' values lying between the 
two curves the homogeneous AF state is self-consistently 
obtained. 




U 



FIG. 6. For negative t' the phase diagram shows no region 
of stability for 8 doped holes per spin (a; « 8%). The AF state 
is stable with respect to longitudinal perturbations only (in 
the region between the two curves). The crossover at 17 ~ 4 
coincides with the vanishing of the AF band gap. 



0.125 
0.120 
0.115 
0.110 
0.105 



0.100 



1 1 1 


1 1 1 1 




N = 14 x 14 




x 8% 


1 1 1 


U = 10 

i i i i 



0.00 -0.05 -0.10 -0.15 -0.20 -0.25 -0.30 -0.35 -0.40 



FIG. 7. For negative t' , the largest eigenvalue A™ ax re- 
mains well above 1/(7, indicating no region of stability for 
the AF state with respect to transverse perturbations. 

However, the response to transverse perturbations shows 
that the AF state is actually unstable. This instability is 
reflected in the behaviour of the largest eigenvalue \™ ax 
of the [x ] matrix, shown in Fig. 7 for the specific case 
of U = 10. It is seen that with increasing |t'|, A™ ax ap- 
proaches 1/U from above, but saturates well above and 
starts increasing again. Similar behaviour is seen for in- 
termediate (U — 5) and weak (U = 3) coupling, and 
also for smaller doping concentration of 4 holes per spin 
(x ;=y 4%). This leads to the conclusion that for negative 
t', the hole-doped t — i'-Hubbard model does not possess 
a stable AF ground state. 



V. TRANSVERSE SPIN FLUCTUATIONS 



The stable mean field AF state for the positive (nega- 
tive) t' model and hole (electron) doping yields a metallic 
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AF state, characterized by intraband (gapless) particle- 
hole excitations, in addition to the interband excita- 
tions. In this section we study the magnon propaga- 
tor, given by x~ + (q,^) = X (<l,u)/1 - U\ ;°(q,w) at 
the RPA level, which represents transverse spin fluc- 
tuations about the broken-symmetry AF state. This 
study should be especially relevant to the electron-doped 
cuprate Nd2- x Ce x Cu04. We focus on the magnon damp- 
ing, which is associated with the imaginary part of 
X°(q, w), and arises from magnon decay into intraband 
particle- hole excitations. Magnon damping is expected 
to play a crucial role in determining the energy depen- 
dence of the local spin susceptibility ^ q Imx~ + (q, w) in 
the doped antiferromagnet, which is probed in neutron- 
scattering experiments. E3] 



A. Evaluation of [x°(q, 



Here we examine the full q, lu dependence of 
[X°(q, w)] intra in the AF state, given in Eq. (8). Con- 
sidering first the imaginary part, we obtain for the AA 
matrix element 

Im[ X °(q, u)YlT = E ^ - ) «W- q " E F ) 



4A 2 



(18) 



Here terms of the type e£/4A 2 are the AF coherence 
factors, and the 9 functions ensure that states k lie be- 
low the Fermi energy while k — q lie above it. The k 
states which contribute in the above expression lie in the 
crescent-shaped region of Fig. 3. 

In the t' » Jx limit considered earlier, and for the 
k states of interest, we have — E^ w e k -q — e k ~ 

At' qa cos 9 and £k « —2t(a 2 /2)cos29 for small q x = q 
and vanishing q y . Substituting in Eq. (18), and taking 
lu > 0, we obtain 



Im[ X °(q, W )]^ a = 27r 



7r/2 ad8 q cos 9 
(2tt) 2 



I A ' V ~2 ^ C ° s2 2 ^ S ^ 4:t ' qa cos ~ w ) 



t 2 Aa 4 
4A 3 8TTt' 



cos 9 
sin 9 



cos 2 29 



(19) 



At' qa cos 9—u 



for lu < At'qa. The delta function never clicks for 
lu > At'qa and the imaginary part vanishes. Thus a new 
energy scale lu x = At ' q a characterizes the imaginary part 
of Im[x (q, u>)], and therefore the spin fluctuation spec- 
trum. 




FIG. 8. Plot of the imaginary part Imx°(q, to) 

vs. u/At'qa, showing the linear dependence on lo for 
lu/ At'qa « 1. 

The lu dependence is particularly simple in the low- 
energy limit lu << lu x . In this limit 9 w tt/2, so that 
sin 9 ps 1 and cos 2 29 m 1 , and therefore we obtain 



Im[ X °(q, W )] 



AA 



t 2 Aa 4 
4A 3 8nt' 



(20) 



which scales like Ax 2 /t' and increases linearly with 



On the other hand, as lu — ► lu x from below, cos 9 — > 1 
in Eq. (19), and therefore Im[x°(q, lu)] diverges like 
At intermediate energies, when lu w 
jj , \ 2. we have 9 w 7r/4 so that the imaginary part 



l/v/T-^ 2 . 

j;/v2, we have 9 
nearly vanishes due to the cos 2 29 term, which originates 
from the AF coherence factor. The full lu dependence is 
shown in Fig. 8. 

We now consider the other matrix elements of 
Im[x°(q, cj)]- For the off-diagonal term we obtain (again 
for lu > 0) 



Im[ X °(q, W )]£fT = J2 ~ ) ^k- q - Ef 

k 

e k £ k-q cfrpQ P 9 , ,\ 



(21) 



The relevant AF coherence factor in this case is 
e k £k- q /4A 2 instead of e£/4 A 2 . Similarly the B- 
sublattice term Im[x°(q, ^)]bb involves the factor 
e£ /4A 2 . Now 6k_ q ~ tk in the limit of small q, there- 
fore to leading order in q all the matrix elements yield the 
same contribution obtained above. Including this intra- 
band contribution to the imaginary part of [x°(q, ^)], wc 
obtain 



[X°(q,^)] 



1 . 

U' 



A 3 



i - i r° + lu 



7q -zr° 
- i r° - lu 



(22) 



Here the energies lu and T° are measured in units of the 
magnon energy scale 2J, and the imaginary term T° is 
given by (i 2 /A 3 )(r n /2 J) = Im[x°(q, w)] intra , given in 
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Eqs. (19) and (20). In the above equation we have re- 
tained only the intra-band contributions to the imaginary 
terms and have neglected those for the real terms, which 
will be small in the limit of small doping concentration. 



B. Magnon Propagator 



Substituting in the RPA expression [x h (q, u>)] = 
[x (q,w)]/l — U[x°(<l, oj)] for the magnon propagator, 
we obtain 



br + (q,^)] = 



(i-*r°-w) -( 7q -<r°) 

-(7q-tr°) (l-iT° + u>) 

1 

> 2 -w 2 -2ir (i- 7q ) 



(23) 



where aj q = W 1 — 7 2 is the bare magnon energy (in units 

of 2 J) for the two-dimensional antiferromagnet. In the 
small momentum limit 7 q ~ 1, so that in the energy 
denominator above the imaginary term 2iT°(l — 7 q ) « 
«r°(l — 7 2 ) = iT°uj^, and the magnon propagator expres- 
sion simplifies to 



br + (q^)] = 



(i-ir°-w) -(7q-<r°) 

-(7q-«r°) (l-tT°+Lj) 



1 



CO ■ 



(i-*r°/2) 

4 (24) 



where 0(r°) 2 terms have been dropped. The magnon 
damping term is finally obtained as T — r°cj q /2. In the 
low-energy limit u << u> x , where T° is given by Eq. (20), 
the magnon damping becomes proportional to energy 



r 

2J 



J'/x 3 / 2 



(25) 



indicating that the characteristic energy scale for magnon 
damping goes as J'/x 3 ^ 2 , where J' — 4i /2 /C7 is the ex- 
change energy corresponding to the NNN hopping t'. 

As discussed earlier, for to > lj x = At'aq, the imagi- 
nary term Tq vanishes, and therefore with an infinitesi- 
mally small damping term (r — > 0) the magnon propa- 
gator in Eq. (24) yields a delta function for the spectral 
function Trlmx '"(q, oj) at lo = w q . If J > At 1 a, then 
cij q > to x for all the small q modes, and in this case the 
local transverse spin susceptibility Trlm% '"(q, lo) ~ 
JqdqS(uj — u>q)/u>q — constant for u> > At' a. How- 
ever, in the low-energy limit {to « to x ), evaluation of 
^ q TrImx h (q, to) from Eq. (24) shows that the local 
spin susceptibility vanishes as uj — > 0. This result is qual- 
itatively similar to the pseudogap behaviour seen in the 
spin excitation spectrum of doped cuprates in neutron- 
scattering experiments, pa] 



VI. CONCLUSIONS 

By significantly modifying the AF band energies near 
the top (bottom) of the lower (upper) band, a positive 
(negative) t' term strongly suppreses the intraband con- 
tribution to the transverse response eigenvalue, resulting 
in a stabilization of the AF state for finite hole (elec- 
tron) doping. We have derived analytical expressions for 
the interband and intraband contributions to the trans- 
verse response eigenvalue in the strong coupling and weak 
doping limit, yielding a cubic equation in t' for the phase 
boundary, solutions of which yield the region of stability 
of the AF state. We find that doping supresses the t'- 
induced frustration due to the competing interaction J' 
and interestingly stabilizes the AF state against the in- 
stability towards the F-AF phase with Q = (vr, 0). These 
analytical results are in quantitative agreement with a 
numerical stability analysis on finite lattices, which also 
yields the full phase diagram in the if — U space, indi- 
cating regions of stability of the AF state with respect to 
both transverse and longitudinal perturbations. 

For negative (positive) t', however, the AF state is 
found to remain unstable for any hole (electron) dop- 
ing. For a fixed sign of t', the location and shape of the 
Fermi surface, the DOS near E-p, and therefore the in- 
traband contribution, are all very different for electron 
and hole doping, as the two AF bands are highly asym- 
metric. This accounts for the strong dependence of the 
stabilization on the type of doping, and qualitatively ex- 
plains why AF ordering in the electron-doped cuprate 
Nd2- x Ce x Cu04 is much more robust than in the hole- 
doped La2- x Sr x Cu04. 

Finally, the study of transverse spin fluctuations in 
the metallic state of the doped AF yields an interesting 
low-energy result T ~ to for magnon damping, arising 
from magnon decay into particle- hole excitations. This 
is of relevance in the context of phcnomcnological theo- 
ries such as the nearly AF Fermi-liquid theory, put for- 
ward to explain the Knight shift experiments and de- 
scribe the non-Korringa temperature dependence of the 
nuclear spin-lattice relaxation rate in doped cuprates. As 
the AF-state energy spectrum is what essentially deter- 
mines the imaginary part of the particle-hole propaga- 
tor, charge fluctuations will also exhibit a similar lu/lu x 
dependence. A detailed study of spin and charge fluctu- 
ations in the metallic AF state, and their consequences 
on the AF long-range order, the local spin susceptibil- 
ity X) q lm X (Qjk-Oj an d the electronic spectrum is cur- 
rently in progress. fl26| 



VII. APPENDIX 

We consider a mean-field AF state with spin polariza- 
tions in the z direction, so that the mean-field Hamilto- 
nian is spin diagonal. Let the eigenvalues and cigenfunc- 
tions in this state be denoted by {£i a , \lcr)}. The state 
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\la) has amplitudes <j>\ a on site i for spin a and zero ampli- 
tude for the opposite spin. Then the first-order correction 
to the state \l J.) due to a small transverse perturbation 
A i is obtained as 



\5l T)=E 

m T 



|mt)(mT |A X |U) 

£ll - £mt 



(26) 



The correction |<52 f) yields spin-f amplitudes generated 
by the perturbation to the state \l I), which has only 
spin-| amplitudes. A similar equation yields the correc- 
tion \Sl i) for the \l f) state. 

Now, the transverse mean-field potential Aj_ = 
— U (at ) generated due to the transverse perturbation 
can be derived from the exact eigenstates \L), and we 
obtain 





1 



(27) 



to first order in the perturbation. Here the ' indicates 
that the sum is over all occupied states. Substituting for 
the amplitudes, we obtain 



-^EE 



Sit — £; 



A' 



mi 



-'ml 



(28) 



As the two terms above are antisymmetric under ex- 
change of labels I and m, all terms with states m below 
the Fermi energy cancel pairwise, and we obtain 




Hence the transverse response matrix [A] = [x°], the bare 
antiparallel-spin particle-hole propagator in the static 
limit (u> = 0). 



9 
10 



11 
12 



18 
19 
20 

21 
22 
23 



Electronic address: avinas@iitk.ac.in 
Present address: Laser Physics Division, Centre for Ad- 
vanced Technology, Indore - 450213, India. 
H. R. Krishnamurthy, C. Jayaprakash, S. Sarker, and W. 
Wenzel, Phys. Rev. Lett. 64, 950 (1990); 
S. Sarker, C. Jayaprakash, H. R. Krishnamurthy, and W. 
Wenzel, Phys. Rev. 43, 8775 (1991). 
Z. Y. Weng and C. S. Ting, Phys. Rev. B 42, 803 (1990). 
A. Singh, Z. Tesanovic, and J. H. Kim, Phys. Rev. B 44, 
7757 (1991). 

A. Singh, Z. Tesanovic, and J. H. Kim, Pramana - J. 
Phys. 38, L211 (1992). 

This AF ground state would be exact for the generalized 

Af orbital Hubbard model in the Af — * oo limit, as all 

fluctuations involving powers of 1/AT would vanish. See, 

e.g., A. Singh, Phys. Rev. B 43, 3617 (1991). 

M. S. Hybertsen, E. B. Stechel, M. Schliiter, and D. R. 

Jennison, Phys. Rev. B 41, 11068 (1990). 

S. Bacci, E. Gagliano, R. Martin, and J. Annett, Phys. 

Rev. B 44, 7504 (1991). 

C. Almason and M. B. Maple, in Chemistry of High- 
Temperature Superconductors, edited by C. N. R. Rao, 
World Scientific, Singapore (1991). 

T. Sasaki et ai, Solid. St. Commun. 75, 93 (1990). 

S. Carter et ai, Phys. Rev. Lett. 67, 3440 (1991); 

W. Bao et ai, ibid. 71, 766 (1993). 

S. Sudo, J. Mag. Mag. Mat. 114, 57 (1992). 

A. J. Millis, H. Monien, and D. Pines, Phys. Rev. B 42, 

167 (1990). 

P. Benard, L. Chen, and A. -M. S. Tremblay, Phys. Rev. 
B 47, 589 (1993); 

A. Veilleux, A. Dare, L. Chen, Y. M. Vilk, and A. -M. S. 

Tremblay, Phys. Rev. B 52, 16255 (1995). 

T. Tohyama and S. Maekawa, Phys. Rev. B 49, 3596 

(1993). 

G. Stemmann, C. Pepin, and M. Lavagna, Phys. Rev. B 
50, 4075 (1994). 

O. K. Andersen, A. I. Liechtenstein, O. Jepsen, and F. 
Paulsen, J. Phys. Chem. Solids 56, 1573 (1995). 
A. V. Chubukov and K. A. Musaelian, J. phys., Condens. 
Matter, 7, 133, (1995). 

D. Duffy and A. Moreo, Phys. Rev. B 52, 15607 (1995). 
D. Duffy and A. Moreo, Phys. Rev. B 55, 676 (1997). 
W. Hofstetter and D. Vollhardt, Ann. Phys. (Leipzig) 7, 
48 (1998). 



A. Singh, Report No. cond-mat /Oil 2442 



A. Singh, Report No. cond-mat /0006079 (2000) 



(2001). 

Condens. Mat- 



H. Ghosh and R. R. dos Santos, J. Phys.: 
ter 11, 4499 (1999), and references therein. 
A. Singh and Z. Tesanovic, Phys. Rev. B 41, 614 (1990); 
Phys. Rev. B 41, 11 457 (1990). 

H. F. Fong, P. Bourges, Y. Sidis, L. P. Regnault, J. Bossy, 
A. Ivanov, D. L. Milius, I. A. Aksay, and B. Keimer, 
Phys. Rev. B 61, 14773 (2000). 
A. Singh, to be published. 



10 



